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Abstract 

Traditional theory of many-electron atoms and ions is based on the 
coefficients of fractional parentage and matrix elements of tensorial 
operators, composed of unit tensors. Then the calculation of spin- 
angular coefficients of radial integrals appearing in the expressions of 
matrix elements of arbitrary physical operators of atomic quantities 
has two main disadvantages: (i) The numerical codes for the calcula- 
tion of spin-angular coefficients are usually very time-consuming; (ii) 
f-shells are often omitted from programs for matrix element calculation 
since the tables for their coefficients of fractional parentage are very 
extensive. 

The authors suppose that a series of difficulties persisting in the 
traditional approach to the calculation of spin-angular parts of matrix 
elements could be avoided by using this secondly quantized methodol- 
ogy, based on angular momentum theory, on the concept of the irre- 
ducible tensorial sets, on a generalized graphical method, on quasispin 
and on the reduced coefficients of fractional parentage. 
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1 Introduction 



Modern atomic spectroscopy studies the structure and properties of prac- 
tically any atom of the periodic table as well as of ions of any ionization 
degree. Particular attention is paid to their energy spectra. For the in- 
vestigations of many-electron atoms and ions, it is of great importance to 
combine experimental and theoretical methods. Nowadays the possibilities 
of theoretical spectroscopy are much enlarged thanks to the wide usage of 
powerful computers. Theoretical methods utilized must be fairly universal 
and must ensure reasonably accurate values of physical quantities studied. 

Many-electron atom usually is considered as many-body problem and is 
described by the wave function constructed from the wave functions of one 
electron, moving in the central nuclear charge field and in the screening field 
of the remaining electrons. Then the wave function of this electron may 
be represented as a product of radial and spin-angular parts. The radial 
part is usually found by solving various modifications of the Hartree-Fock 
equations and can be represented in a numerical or analytical forms (Froese 
Fischer 1977) whereas the angular part is expressed in terms of spherical 
functions. Then the wave function of the whole atom can be constructed in 
some standard way (Cowan 1981, Jucys and Savukynas 1973, Nikitin and 
Rudzikas 1983) starting with these one-electron functions and may be used 
further on for the calculations of any matrix elements representing physical 
quantities. 

During the last two decades a number of new versions of the technique 
(so-called Racah algebra) to cope with spin-angular parts of the wave func- 
tions and matrix elements have been suggested (Rudzikas 1991). Among 
them the second quantization and quasispin techniques turned out to be 
of particular efficiency (Judd 1967, Rudzikas and Kaniauskas 1984). The 
usage of graphical methods (Jucys and Bandzaitis 1977) allowed one to find 
general expressions even for rather complex cases of matrix elements. All 
this enabled one to formulate fairly consistent and general non-relativistic 
and relativistic theory of many-electron atom and processes of its interac- 
tion with electromagnetic radiation (Rudzikas 1996). The abovementioned 
methods are applicable for the both variational and perturbative approaches 
for various coupling schemes of spin and orbital momenta. 

Practically we have to solve so-called eigenvalue problem 

HV = EV, (1) 
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where \I/ is the wave function of the system under investigation and H is 
its Hamiltonian. In various versions of perturbation theory such equation 
usually serves as the starting point for the further refinements. It turned 
out that for very large variety of atoms and their ionization degrees the 
so-called Hartree-Fock-Pauli Hamiltonian leads to highly accurate energy 
values (Nikitin and Rudzikas 1983, Rudzikas 1996) that is why it is widely 
used in many methods and computer codes. 

In order to calculate the energy spectrum of atom or ion we have to find 
the expressions for the matrix elements of all terms of the Hamiltonian con- 
sidered. For complex electronic configurations, having several open shells, 
this is a task very far from the trivial one. For the optimization of their ex- 
pressions one has to combine the methods of the angular momentum theory, 
irreducible tensorial sets, tensorial products in a coupled form, coefficients 
of fractional parentage with the utilization of the graphical (diagrammatic) 
methods, second quantization and with the accounting for the symmetry 
properties of the system under consideration in the additional spaces, for 
example, quasispin space. This paper describes one such possibility. 

Unfortunately, practical calculations show that all realistic atomic Hamil- 
tonians do not lead straightforwardly to eigenvalue problem (1). Actually we 
have to calculate all non-zero matrix elements of the Hamiltonian considered 
including those non-diagonal with respect to electronic configurations, then 
to form energy matrix, to diagonalize it, obtaining in this way the values 
of the energy levels as well as the eigenfunctions (the wave functions in the 
intermediate coupling scheme). The latter may be used then to calculate 
electronic transitions as well as the other properties and processes. Such 
necessity raises special requirements to the theory. 

The total matrix element of each term of the energy operator in the 
case of complex electronic configuration will consist of matrix elements, de- 
scribing the interaction inside each shell (in relativistic case - each subshell) 
of equivalent electrons as well as between these shells. Going beyond the 
single-configuration approximation we have to be able to take into account in 
the same way non-diagonal, with respect to configurations, matrix elements. 
Starting at the very beginning with the second quantization and quasispin 
methods we are in a position to fulfill all these requirements. Below we shall 
describe the approach suggested in the more details. 
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2 Tensorial Form of the Operators 

According to the method of second quantization (Judd 1967, Rudzikas and 
Kaniauskas 1984) any one-particle operator 

F = J2a l a+(i\f\j), (2) 
can be expressed in the following tensorial form: 

(Kcr)r 



x a 



F= Yl F(hj)= E [«,cr]- 1/2 (n,A i ||/^)||n J A j 

(3) 

where i = nilism^nis^ A = Is, [k, a] = (2k + 1) (2a + 1) , (riiXiWf^^WrijXj^J 
is the one-electron submatrix (reduced matrix) element of operator F, and 
is the electron creation operator. The tensor a 3 is defined as 

where a^it ^ s the electron annihilation operator. From tensorial point of 

view it is better to consider tensor a as electron annihilation operator 

(Ko-)r 

(see Section 4) . The product of tensors 



a^x a 



denotes tensorial 



part of operator F. Here the rank k of orbital space is coupled to the spin 
space rank a to form a tensorial product of total spin-angular rank T. As we 
shall see, this expression is very effective for the calculation of spin-angular 
coefficients for any one-particle operator. This expression is a general one 
and the tensorial form of any one-particle physical operator may be obtained 
from it. For example, the spin-orbit interaction operator has the tensorial 
structure re = l,cr = l,r = 0, and its submatrix element is 

fx x 1 / 2 
(^AII/i-JlM,) = ~ za2 & + !) (2*< + 1)J (n^jl/r 3 !™,/,) Sfalj) 

' (5) 

Any two-particle tensorial operator 

G = ^ E a i a j a f a i> (hj\9\i',f) (6) 
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can be expressed in two well-known forms (Rudzikas and Kaniauskas 1984). 
In the first form the operators of second quantization follow in the normal 
order: 



Gi 



E E E(- 1 ) fe_P [^2,K' l2 ,a 12 ,a[ 2 ] 1/2 x 

iji i i 



x (n i \ i n j X j \\g^ K2k ^ a ^\\n' i X' i n' j X' j 



x < 



Ki K 2 k 



> < 



S S cr 12 
(7l (7 2 k 
S S <J\ 2 



> x 



(«12CTl 2 ) 



x a 



( K 'l2 CT 12) 



(kk) 



p-p 



(7) 



where (njAjnjAj||g( K1K2fc ' CT1CT2fc )||n^A^-A^ is the two-electron submatrix el- 
ement of operator G. 

In another form the second quantization operators are coupled by pairs 
consisting of electron creation and annihilation operators. In tensorial form: 



G u 



E G H (iji'j') 



ndi,ndj,n f .V,n f V. 



E E (-!) fc " P (niA i n i A i ||^*'^*)||^n;A;) 



x{[K l7 K 2 ,ai,a 2 ] 1/2 



,(A0 



x a 



(a')1 ( KlfJ1 ) 



a v '^x a 



(a'.)1 ( K2 °" 2 ) 



(-1) 



ii+i' J K l K 2 k \ J CJi (7 2 & 



/' /• /• 

fc^ fcj 



a^x a 



(*;) 



(fcfc) 



p-p 



(8) 

The expression (7) consists of only one tensorial product whereas (8) 
has two, but the summation in the first formula is also over intermediate 
ranks k± 2 , &12, k' 12 and a' 12 , complicating in this way the calculations. The 
advantages or disadvantages of these alternative forms of arbitrary two- 
electron operator may be revealed in practical applications. 
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In these forms the product of second quantization operators denotes ten- 
sorial part of operator G. For instance, the tensorial structure of electrostatic 
(Coulomb) electron interaction operator is the same as that of orbit-orbit 
interaction, K\ = K2 = k, o\ = 02 = (Jucys and Savukynas 1984), and only 
the two-electron submatrix elements ^njAjnjAj||g( KlK2fc ' CT1CT2fc )||n^njA^ of 
these operators are different. In the case of electrostatic interaction: 

(niXiUjXj I \g { couiomb I WiKn'j = 
= 2 [kf 2 (li\\CW\\l^ (IjWC^WI'j) R k (nikn^njljn'/j) . (9) 

From (9), by (7) and (8), we finally obtain the following two secondly 
quantized expressions for Coulomb operator: 



Vi 



1 y ^_^yj+ii+k+K 12 



[KIS,^] 1 / 2 (li\\CW\\ty X 



(K12IT12) 



x a 



(K120"12)' 



(00) 



E , E {hWC^Wl'i) (^||C( fe )||/;.) R k (nMl'i.rijljn'/j) 

nilin-ihn'l'.n f .l f . k 



(10) 

x 



a^x a 



x{[fc]- 1/2 



( A 



(fcO) 



a^x a 



x 



a ( A i)x a 



(koy 



(00) 



(00) 



+ 



Sinjl^n'^)}, 



(11) 

The tensorial expressions for orbit-orbit and other physical operators in 
second quantization form may be obtained in the same manner. 

It is worth mentioning that the expressions (10) and (11) embrace, al- 
ready in an operator form, the interaction terms both the diagonal ones, 
relative to configurations, and the non-diagonal ones. Non-diagonal terms 
define the interaction between all the possible electron distributions over the 
configurations considered, differing by quantum numbers not more than two 
electrons. 
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The merits of representing operators in one form or another (10) or (11) 
are mostly determined by the technique used to find their matrix elements 
and quantities in terms of which they are expressed. 



3 Generalized Graphical Method 

In this paragraph we shall sketch the generalized version of graphical tech- 
nique, in which not only one- and two-particle operators are presented in 
tensorial form (such graphs are analogical to Feynman-Goldstone diagrams 
but they do not depend on magnetic quantum numbers (Merkelis et al 1986a, 
b), but which allows also to represent graphically any tensorial product of 
the second quantization operators and to perform graphically the opera- 
tions with the secondly quantized operators as well as with their tensorial 
products (Gaigalas et al 1985, Gaigalas 1985, Gaigalas and Merkelis 1987). 
Such graphical technique is most suitable to represent any one- and two- 
particle operator already presented in tensorial form and to found general 
expressions for their matrix elements. 

In this methodology the item under summation sign of the one-particle 
operator (3) has the following graphical form: 



F(i,j) = Ai = [K,a]~ l/2 (niXiWf^WnjXj) A 2 A 3 , 



(12) 



where the diagrams A±, Ai and A3 are presented on Figure 1. As we see, the 
diagram of the operator itself, namely A±, is similar to the usual Feynman- 
Goldstone diagram (Lindgren and Morrison 1982), although here the sum- 
mation over magnetic quantum numbers m\ is performed. The product of 
the diagrams A2, A3 represents the tensorial structure of the operator: 



A 2 A 3 = 


a( A «)x« (y ' 


(Kcr)r 

= E 




(kct) 






mr m K ,m a 




m K m a 



k a r 

m K m a m r 



where A2 equals: 



Ao = 



(k<t) 



(13) 
(14) 



whereas A3, by (Jucys and Bandzaitis 1977), is equal to: 
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Figure 1: Diagrams for one-particle operators. 
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A, 



K 

m K 



rriY 



(15) 



The heavy line in the diagram A3 represents the resultant momentum 
r whereas the plus sign of the vertex means that the momenta k and a 
are coupled into the resultant T in counter-clockwise direction. From the 
symmetry properties of the Clebsch-Gordan coefficients the equality follows: 



A 3 = (-1) K+ ^ T A 4 . 



(16) 



Then we can conclude that, if we change the sign of any vertex, then the 
phase multiplier of the form ( — l) K+<T_r occurs. 

The electron creation operator a^*) has the following graphical form 
(Figure 1, A 5 ): 



(17) 



whereas a 



Ar. 



(18) 



Thus, it is obvious that the diagram A2 consists of the second quantiza- 
tion operators a { i] and a as well as of the Clebsch-Gordan coefficients 



A 7 



u 



h 3 



K 

m 



K 



A* 



s 

m s 



s 

m s 



a 



(19) 



which couple these operators into tensorial product and which may be ob- 
tained from the diagram A2 if to omit in them the graphical symbols of 
the second quantization operators. It is necessary to bear in mind that, 
while writing down the algebraic expression from the diagram A 2 , always in 
tensorial product there must be in the first place the second quantization 
operator, which is above the vertex "a", whereas the second place must oc- 
cupy the operator, which is below the vertex "a" in the diagram Ai- The 
scheme of their coupling into tensorial product is defined by the sign of the 
vertex. 

The first form (7) of two-particle operator Gi(iji'j') is represented by 
the following diagram (Figure 2, B\): 
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kk 



rijXj 



n>\> 3 



kk 



rijXj 



n i X l 
+ 



n'jX'j rijXj 



kk 



n'X 



Bo 



B, 




B4 -B5 Bq 



Figure 2: Diagrams for two-particle operators. 
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G I (iji'j') = B 1 =-\ £ 2 (-l)*"^^,^, (71, (7 2 
Kl2« 1 2« T 12CT 12 



1/2 



3 




K 'l2 




^2 


jfe 






«12 



> x 



x < 



s s a 12 

(7i CJ 2 /C 
S S CJ12 



>B 3 



(20) 



whereas the second (8): 



G n (iji'f) = B 2 + B A = 
p 



(21) 



We emphasize here that the winding line of interaction in the Feynman- 
Goldstone diagram corresponds to the operators of second quantization in 
the normal order (Figure 2, B\). Whereas the dotted interaction line indi- 
cates that the second quantization operators are ordered in pairs of creation- 
annihilation. In the latter case first comes the pair on the left side of a 
Feynman-Goldstone diagram (Figure 2, B2). Such a notation of two kinds 
for an interaction line is meaningful only for two-particle (or more) opera- 
tors, since for any one particle operator both the winding and dotted lines 
correspond to the same order of creation and annihilation operators. 

From expressions (20), (21) we see that the two-particle operator in the 
first form is represented by one Feynman-Goldstone diagram B±, whereas in 
the second - by two diagrams B2 and B4. The diagrams, corresponding to 
tensorial product, have the following algebraic expressions: 



5, 



a (Ai) x a (\j) 



(K12CT12) 



~(AJ) ~(A') 



( K 'l2°"l2) 



a (Al) x a" 



(a')1 ( Kicri ) 



a (Aj) x a' 



(a')1 ( K2fT2 ) 



(kk) 

P~P 
(kk) 

P~P 



(22) 
(23) 
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o (Ai) x a 



(fcfc) 



(24) 



Thus, the obtaining of algebraic expressions from the diagrams B%, B§ 
and Bq is similar to the case of the diagram A<i- The positions of the 
second quantization operators in the diagram define their order in tensorial 
product: the first place in tensorial product occupies the upper right second 
quantization operator, the second - lower right, after them the upper left 
and lower left operators follow. The angular momenta diagram defines their 
coupling scheme into tensorial product. 

Thus, obeying these rules it is possible to easily find the algebraic coun- 
terparts of the diagrams, not forgetting that the arrangement of the oper- 
ators must not contradict to their coupling order, i.e. only neighbouring 
second quantization operators are coupled into tensorial product and their 
disposition order corresponds to coupling scheme. Otherwise some graphical 
operations are necessary. Let us present the simplest of them below as the 
example for the case, when we have to change the disposition of the second 
quantization operators and coupling scheme in the tensorial product. 

Suppose, we have the following correspondence between diagrams (Fig- 
ure 3): 

Ci — C 2 , (25) 
in which the second quantization operators are in the order 

a (A 3 ) q(A 4 ) fl (Ai) 

a^ 2 ^. Our goal is to obtain the diagram corresponding to the order a^ Xl ^ 
q(A2) a (A 3 ) q(A4)_ Bearing in mind that the second quantization operators 
anticommute with each other and they all act on different electronic shells 
and we are not changing the order of their coupling into tensorial product, 
we arrive at 



d — (-1) 4 C 3 = C 3 . 



(26) 



Let us also discuss another situation: we have defined the disposition 
of the operators and we want to couple them into certain tensorial prod- 
uct. Suppose that we want to represent graphically the following tensorial 
product: 



x a 



(A 2 )] ^ x 



Q (As) x g(A4) 



(K2C2) 



(kct) 



(27) 



12 



ra 3 A 3 



kk 



n 2 A 2 



n 3 A 3 
+ 



niAi 



n 2 A 2 



Ci 



Co 



niAi 
+ 
n 2 A 2 



+ 



K\ (T\ K2O2 

kk 



nzX'i 



714 A4 



c 3 



+ 
n 2 A 2 



— r — 

kk 



n3X 3 



71^4 



Figure 3: Diagrams for graphical transformations. 
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For this purpose we have to rearrange the generalized Clebsch-Gordan 
coefficient, which is defining the scheme of coupling of the operators into the 
tensorial product. It is easy to notice that this coefficient will properly define 
the tensorial product, if we change the sign of the vertex "a" in diagram C3. 
Making use of (16) we find: 

Ci — ► (_i)«i+«2-«+ CT i+ CT2 - CTC < 4 _ (28) 

The procedures described are fairly simple, however, they are sufficient 
for the majority of cases. The more complete description of this generalized 
graphical approach may be found in Gaigalas et al 1985, Gaigalas 1985, 
Gaigalas and Merkelis 1987. 



4 Quasispin Formalism 

A wave function with u shells in LS coupling may be denoted in the form 



?p u (LSM l M s ) = 
= \n 1 l? 1 n 2 l^...n u l^a 1 L 1 S 1 a 2 L 2 S 2 ...a u L u S u ALSM L M s ), 



where A stands for all intermediate quantum numbers, depending on the 
order of coupling of momenta ajLjS'j. 

As we shall see later on, it is very convenient for the calculations of matrix 

elements to use quasispin formalism. Then am\ and a are components of 

the tensor amqm x j having in additional quasispin space the rank q = ^ and 

projections m q = ±\, i.e. af^ = a^}m s and 0^1^ = a L;L s - 

In the quasispin representation, for a wave function of the shell of equiv- 
alent electrons \nl N aLS) a label Q - quasispin momentum of the shell - is 
introduced, which is related to the seniority quantum number v, namely, 
Q = (21 + 1-v) /2, and its projection, M Q = (N — 21 — 1) /2. Here a 
denotes all additional quantum numbers needed for the one-to-one classifi- 
cation of the energy levels. Then, we can rewrite the wave function (29) 

^ u (LSM L M S ) = 
= \n 1 l 1 n 2 l 2 ...n u l u a 1 L 1 SiQ 1 M Ql a 2 L 2 S 2 Q 2 MQ 2 ... 

a u L u S u Q u M Qu ALSM L M s ). (30) 

Making use of the Wigner-Eckart theorem in quasispin space of a shell 

l N , 
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(l aQLS\\\T^\\\l a'Q'L'S') (31) 



(l aQLSM Q \\T^ s) \\l a'Q'L'S'M' Q ) = (-l) 2q [Q]' 1 ' 2 x 

Q' q Q 

M' Q m q Mq 

it is possible to define the notions of a completely reduced matrix element 
(l aQLS\\\T^\\\l a'Q'L'S') and subcoefficient of fractional parentage (re- 
duced coefficient of fractional parentage) (l aQLS\\\a^ qls ^\\\l a'Q'L'S'). In 

(31) is any tensor with rank q and its projection m q in quasispin space 

and on the right-hand side of this equation only the Clebsch-Gordan coeffi- 
Q' q Q 



cient 



M Q 



M Q 



depends on the number N of equivalent electrons. 



According to Rudzikas and Kaniauskas 1984 we have the following re- 
lation between the coefficients of fractional parentage and completely re- 
duced matrix elements (l aQLS\\\a^\\\l a'Q'L'S') of the operator of sec- 
ond quantization a 



(qls). 



(l N aQLSW^' 1 (a'Q'L'S')l) 



Q' 1/2 Q 
M' Q 1/2 M Q 



(l aQLS\\\a^\\\l a'Q'L'S') . (32) 



Tables of numerical values of (l aQLS|||a(9 ia )|||i aQ'L'S') are presented 
in Rudzikas and Kaniauskas 1984 when I = 0, 1, 2. For the tensorial product 
of two one-electron operators, the submatrix element equals 



nl N aQLS\\ 

E {Qr 1/2 



(gX) 

0-m q i x a m q 2 



(kik 2 ) 



\\nl N ' a'Q'L'S' 



e,m e 



q q e 
m q i m q2 m e 



Q' 



Q 

Mq 



x (nl aQLS\\\W( eklk ^\\\nl a'Q'L'S') 



(33) 



On the right-hand side of equations (32) and (33) only the Clebsch- 
Q' e Q 



Gordan coefficient 



MA 



Mq 



depends on the number N of equiva- 



lent electrons. 

(nl aQLS\\\W( eklk2 )\\\nl a'Q'L'S') denotes reduced in quasispm space 
submatrix element (completely reduced matrix element) of the triple ten- 
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sor W^ eklk2 ^ (nl,nl) = a^ 1 ^ x a^ 1 ^ ^ 1 2 \ It is related to the completely 
reduced coefficients (subcoefffcients) of fractional parentage in a following 
way: 

(nl aQLS\\\W^ k ^\\\nl u'Q'L'S') = 
= { _ l)Q+ L + s+Q>+L> + s> + e +kl+ k 2 [eMM] V> x 

x ^ (l aQLS\\ \a (qls) \\\l a"Q"L"S") x 

a"Q"L"S" 

x (l a"Q"L"S"\\\a^\\\l a'Q'L'S 1 ) x 

X { Q' Q Q" } { L> L L" } { 5' S S" } ' (34) 

So, by applying the quasispin method for calculating the matrix elements 
of any operator, we can use the reduced coefficients of fractional parentage 
or the tensors (for example W^ eklk2 ^ (nl,nl)), which are independent of the 
occupation number of the shell for a given v. The main advantage of this 
approach is that the standard data tables in such a case will be much smaller 
in comparison with tables of the usual coefficients and, therefore, many 
summations will be less time-consuming. Also one can see that in such an 
approach the submatrix elements of standard tensors and subcoefficients 
of fractional parentage actually can be treated in a uniform way as they 
all are the completely reduced matrix elements of the second quantization 
operators. Hence, all methodology of calculation of matrix elements will be 
much more universal in comparison with the traditional one. 



5 Matrix Elements in the Case of Two Open Shells 
of Equivalent Electrons 

The aim of this section is to illustrate the usage of abovementioned method- 
ology to obtain the expressions for matrix elements of a two-particle opera- 
tor, when the wave function (30) has two open shells niZf 1 n 2 Z^ 2 . Then it 
may be written as 



V> 2 {LSM L M S ) = \n l l l n2l2a 1 L 1 S 1 Q l MQ 1 a 2 L 2 S 2 Q2M Q2 LSM L Ms). (35) 

To find numerical value of physical quantity of two-electron operator 
one ought to have the expressions for its matrix elements within each shell 
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of equivalent electrons and between each pair of the shells, including non- 
diagonal, with respect to configurations, matrix elements. 

While calculating the diagonal matrix elements between functions (35), 
the quantum numbers rijAj, n'^X 1 ^ rijXj, n'jXj in two alternative expressions 
(7) (8) acquire the following values: 

1. riiXi = n! i X\ = rijXj = n'jXj = n,\l\s. (All the operators of second 
quantization act upon the first shell). 

2. riiXi = n! i X\ = njXj = n'jXj = n 2 l 2 s. (All the operators of second 
quantization act upon the second shell). 

3. riiXi = n'^ = nilis, njXj = n'jX'j = n 2 l 2 s. 

4. rijXj = n'jXj = n\l\S, riiXi = n'^ = n 2 l 2 s. 

5. riiXi = n'jX'j = nil±s, n^X^ = rijXj = n 2 l 2 s. 

6. n^X',1 = njXj = nil±s, rijAj = n'jXj = n 2 l 2 s. 

In the first case the matrix elements of operator in the first (using (7)) 
and the second (using (8)) forms are equal respectively 



{n 1 l 1 ? 1 n 2 l% 2 a 1 L 1 S 1 Q 1 M Ql a2L 2 S 2 Q 2 M Q2 LSM L Ms\G I (lin) 
|n 1 /f 1 n 2 / 2 V2 «' 1 L' 1 51Q' 1 MQ 1 a' 2 L' 2 5 2 Q 2 MQ 2 L'^MiM^) = 

= ~\ E £(-l) fc - p [Ki,K2,ai,ai I ] 1/2 x 



K12K 12 a-i2(T 12 




r (1,1 i\, \\ (^120-12) r~(Ai) ~(Ai) 

a (M) x a (M) x a x a 




n 1 l^n 2 l^a[L[S[Q[M Ql a' 2 L' 2 S 2 Q' 2 M Q2 L'S'M 



-ip-p 



(36) 
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{n l l^n2l2 2 aiL 1 S 1 QiM Ql a2L 2 S2Q2MQ 2 LSM L M s \G n {llll) 
|n 1 Zf 1 n 2 ^ 2 a' 1 L' 1 ^Q' 1 M Ql a' 2 L' 2 ^Q' 2 MQ 2 L^'M[M^) = 

= \Y,i- l t~ P (niAiniAill^^^^^HniAiniAi) x 



x{[ki,K 2 ,(7i,(72] 



'1/2 



:(n 1 l? 1 n 2 l^a 1 L 1 S 1 Q 1 M Ql a2L 2 S2Q2M Q2 LSM L M s \ 



x a 



.(AO 



(fClCTl) 



X 



(Ai) 



x a 



.(Ai) 



(K2CT2) 



(fcfc) 



P-P 



|nilf 1 n2Z^y i L' 11 SiQ / 1 Mg 1 a , 2 ^^Q , 2 Mg 2 L' 1 S , M{ ; M^)- 

Kl K 2 k 1 f <Tl CT 2 fe 1 x 
^1 ^1 ^1 J 1 « S S J 

x{n 1 l^n2l^a 1 L 1 S l Q 1 M Ql a 2 L2S2Q2M Q2 LSM L M s \ 



AM) 



(kk) 



P-P 



|n 1 /f 1 n 2 / 2 W2 aiLi^QiMQ 1 a 2 L 2 5 2 Q 2 MQ 2 L'5'M[M^)}. 
Schematically these expressions may be written down as 



(37) 



{n l l^n2^ 2 a l L l S l Q l M Ql a2L2S2Q2M Q2 LSM L M s \G{llll) 
\n 1 l^n2l^ a [L[S[Q[M Ql a'2L' 2 S' 2 Q'2M Q2 L'S'M' L M' s ) = 
= ^2 6 (Ki 2 ,cJi 2 ,K' 12 ,cJi 2 ,ni, Ai) x 



K>12,&12,k' 12 <7' 12 



x(n 1 l? 1 n2l2 2 a 1 L 1 S 1 Q 1 M Ql a2L2S2Q 2 M Q2 LSM L M s \ 

A^-l (Ki2,cri2,Ki 2 ,cri 2 ,ni, Ai) 
|n 1 Z] Vl n 2 Z 2 Y2 a' 1 L' 11 SiQ / 1 M Ql a 2 L 2 S 2 Q 2 M Q2 L' 1 S / M[M^), 



(38) 



where (ki 2 , <ti 2 , k/ 12 , a^ 2 , ni, Ai) is proportional to the radial part of an 
operator, and A p kk p («i 2 , o"i 2 , k' 12 , cr 12 , ni, Ai), - to the spin-angular part of 



it. In the first form 



Ap k _l {K 12 , £712 , k'i 2 » °12 , n l , Al ) 



a (Ai) x a (Ai) 



(K12C12) 



~(Ai) ~(Ai) 

a x a 



( K 'l2 <T 12) 



(fcfc) 



p-p 



(39) 



whereas in the second form (ki 2 = m, a±2 = ai, k' 12 = 1^2, £7^ 2 = a 2 ) 
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A p kk) (Ki,<7i,K2,02,ni,Ai) = {[Ki,K 2 ,a 1 ,a 2 \ 1/2 x 



x a 



-(Ai)' 



(/tlCTl) 



x a 



^(Ai) 



(kk) 



,(Aii) 



x a 



-(Ai) 



P-P 
-i (fcfc) 



}■ 



(40) 



I ki k 2 k M cri (j 2 /c 
1 /i Zi Zi | I s s s 

So, in order to calculate the spin-angular parts of matrix elements of 
operators (7), (8), we have to obtain at first the matrix elements of operators 

A p kk p (K12, 012, k' 12) Oi2> n i> By us i n g the Wigner-Eckart theorem, we 
find: 

(n 1 l^n 2 l^a 1 L 1 S 1 Q 1 M Ql a 2 L 2 S 2 Q 2 M Q2 LSM L M s \ 

A p kk p (Kl2, 012, K 'l2> °12, n l> A l) 

|n 1 /f 1 n 2 Z 2 V2 a' 1 L' 1< S[Q' 1 M Ql a 2 L 2< S 2 Q 2 MQ 2 L'5'M£M^) = 
L' k L S' k S 

My p M L M S : -p M s 
x(n 1 l^n 2 l% 2 a 1 L 1 S l Q 1 M Ql a 2 L 2 S 2 Q 2 M Q2 LS 

A {kk) («12 , 012 , «'l2 » °12 . n l > A l) 

||n 1 Zf 1 n 2 Z 2 V2 a' 1 L' 1< S[Q' 1 M Ql a 2 L 2< S 2 Q 2 MQ 2 L' < S'). 
Then we proceed with analyzing the submatrix elements. As the oper- 
ator J 4( fcfc ) (k 12 , £7i2, «i2> °12> n i> A ac t s here upon the first shell only, then, 
using the expression (4.7) from Jucys and Savukynas 1973, namely, 



= [L,S]- 1 / 2 



(41) 



{ai3ia 2 j 2 j\\A { i ] \\aijia 2 j 2 j) = 5 (a 2 j 2 ,a' 2 j 2 ) ( 



we obtain: 



x[j,if /2 (amPniaiii) 



(*)i 



Ji 




jl+j2+j'+fe 



(niZ^naZ^aiLiSiQiMg^sLsSsQsMQ.LSH 
>l( fcfc ) (k 12 , 012, Ki2, &12i n l) Al) 
| \n 1 l^n2l2 W2 a[L f 1 S f 1 Q' 1 M Ql a f 2 L f 2 S 2 Q' 2 M Q2 L f S f ) = 

= (_i^L 1 +S 1 +L 2 +S 2 +L>+S'+2k y L g jj ^/j 1/2 x 

J Li L L 2 1 f Si 5 ' W l x 

I L' Li fc J \ 5' 5( A: J 

x(ni/f 1 aiQiLiSip( fcfc ) («i 2 , 012, k'i 2 > °"r2, n i> A i 
HniZ^a'iQ'iL^i). 



(42) 



(43) 
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Then there remains only to obtain the formulas for the following subma- 
trix elements: 



(nl N aQLS\\ 



i (A) x a 



(A) 



(kk) 



\\nl N a'Q'L'S') 



(44) 



(nl N aQLS\\ 



(x) 



- (ki<ti) 



X 



(\) ~( A ) 



(K2V2) 



(kk) 



\\nl N a'Q'L'S'), 
(45) 



(nl N aQLS\\ 



(«i20-i 2 ) 



~(A) ~(A) 
a x a 



( K 'l2°12) 



||nP a'Q'L'S'). 



(kk) 



(46) 



Basing ourselves upon the expressions (33), (34), we straightforwardly 
find the value of a submatrix element (44) . The values of submatrix elements 
(45), (46) follow directly from the expression (2.28) of Jucys and Savukynas 
1973 



(aj\\ [A^ x #( fc 2)] (fc) \\ a 'j') = (-iy+i' +k [k] 1/2 x 
E(^ll^ l) H«V)(aVll^)||aY){ *J ) k 

~.ll All I J J J 



(47) 



So we have: 



(nl N aQLS\\ 



a {gX) x a iqX) 



(fticri) 



a m q z x a m(j4 



(K2CT2) 



\\nl N ' a'Q'L'S') = (-tf+S+L'+SW [fe] 
L' L L" ||S" S S" 



(kk) 



a"Q"L"S" 

x(n/ w aQLS|| 
x(nl N " a"Q"L"S" 



t}] (Kiai) \\nl N " a"Q"L"S" 



&m q i x a mq 2 \\nl 

(qX) (,A)1 ("a"*) 

Um q 3 * u 'm q 4 



Ox 



Hn/^' a'Q'L'S') 



(48) 



Schematically we can express the matrix element in the second case, 
when the operators of second quantization act upon the second shell, as 
follows: 
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{n l l^ 1 n2l2 2 a l L 1 S 1 QiMQ 1 a2L 2 S2Q2MQ 2 LSM L Ms\G{2222) 



|n 1 /f 1 n 2 /fa' 1 L' 1 ^Q' 1 MQ 1 4L^^Q^MQ 2 L'^MiM^) = 



X! ® (Kl2,0 - 12,«'l2)°' / 12) n 2, A 2 ) 



X 



K12,C12,Ki2 ,T 12 

x(n 1 /f 1 n 2 / 2 V2 aiLi < SiQiMQ 1 Q2L 2( S 2 Q2M Q2 L5M L M5| 



|fiiZj Vl n 2 Zf a a' 1 L' 1 5ig , 1 M Ql a' 2 L' 2 5^Q' 2 MQ a L' 1 S:'MiM^), 



^Ip^p («12 , CT12 , «12 > °12 ) n 2 > ^2 



(49) 



and find its value by using Wigner-Eckart theorem, expressions (4.9), (2.28) 
from Jucys and Savukynas 1973 as well as (33) and (34). 

Differently from the first and the second cases, in the third (rijAj = 
n^X^ = n\l\8, rijXj = n'jX'j = n 2 / 2 s) and the fourth (rijAj = n^X^ = 
n 2 £ 2 s, rijXj = n'jX'j = n\l\s) cases the first tensorial form (7) is not conve- 
nient for calculating the matrix elements. This is related to the fact that the 
spin-angular part of matrix elements do not have shape of any expression 
below: 

(n 1 l^ 1 n2l2 W2 a 1 L 1 S 1 Q 1 MQ 1 a 2 L 2 S 2 Q2MQ 2 LS\\ 



\\n 1 l^n2l^a , 1 L' 1 S l 1 Q l 1 M Ql a , 2L' 2 S , 2Q , 2M Q2 L'S l ), 



A( kk) («12 , <7i2 , K12 » °12 . n l » A l ) 



(50) 



(n 1 l? 1 n2l2 N2 a 1 L 1 S 1 Q 1 M Ql a 2 L2S2Q2M Q2 LS\\ 

A {kk) (Kl2,CJl2,Ki 2 ,Cji 2 ,n 2 , A 2 ) 

||n 1 /f 1 n 2 / 2 V2 a' 1 L' 1< S[Q' 1 M Ql a 2 L 2< S 2 Q 2 MQ 2 L' < S'), 
{n 1 l? 1 n2l2 Y2 a 1 L 1 S 1 Q 1 M Ql a 2 L2S2Q 2 M Q2 LS\\ 

' A (K 12 a 12 ) ( niAi ) x B {n' X2 c' X2 ) ( n2 A 2 )] (fcfc) 

||nii^n 2 i^aiA5iQiM^ 1 a' 2 L' 2 5^Q' 2 M^ 2 L'5'), 



(51) 



(52) 



(niZfwf^iSiQiMQ^^^Mc^LSII 

" j4 («12-12) ( n2A2 ) x jB («' 1 2- , 12 ) (mAl)]^ 

||n 1 /f { n 2 ^a / 1 L' 1< SiQ / 1 M^ i a 2 L 2< S 2 Q 2 M^L' < S / ). 



(53) 



Here A^ 12 " 12 " 1 (nX) and ijO'ia' 7 ^) (nA) represent any tensorial operator. 
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Only these shapes (50), (51), (52), (53), in the case of two open shells, 
guarantee the effective use of Racah algebra. That includes the determina- 
tion of zero matrix elements from triangular conditions (for example in (52) 
5 {Li,L[, k 12 ), 5 (Si, S[, <7i2 ), 5 (L 2 , L' 2 , k' 12 ), S (S 2 , S' 2 , a' 12 )) without explicit 
calculation, the use of tables of standard quantities, and the use of quasispin 
(see Section 4) at last. 

Meanwhile the second form (8) allows one to exploit the Racah algebra 
to its full extent, as the matrix elements for third case 



(nil^n 2 l^aiLiSiQ 1 M Ql a2L 2 S 2 Q 2 M Q2 LSM L Ms\G II (12l2) 
|n 1 /{ Vl n 2 ^ 2 a' 1 L' 1 ^Q' 1 MQ 1 a 2 L 2 5 2 Q 2 M Q2 L'^MiM^ = 
= ^(-l)^[ Kl , K2 ,a 1 ,a 2 ]- 1 / 2 x 
v 

x (niAin 2 A2||£r (KlK2fe ' CTlCT2fc) ||niAin2A 2 ) x 
x(nil^n 2 l^aiLiSiQiM Ql a 2 L 2 S 2 Q 2 M Q2 LSM L M s \ 



x a 



AM) 



(KlO-j) 



a (A 2 ) 



x a 



(«2^2)' 



(kk) 



P~P 



|n 1 Zf 1 n 2 Z 2 V2 a' 1 L' 11 S[Q' 1 M Ql a 2 L 21 S 2 Q 2 M Q2 L' 1 S'MiM^). 
and fourth case 



(54) 



(nil^n 2 l^a 1 LiSiQ 1 M Ql a 2 L 2 S 2 Q 2 M Q2 LSM L M s \G II (2121) 
|niij Vl n2^V 1 Li5igiMo 1 Q , 2 L , 2 ^g , 2 M 02 L , 5 , MiM^) = 

= £(-i) 



\k-p r 1-1/2 w 

1 P kl,K 2 ,«7l,0-2 ' X 



x (n 2 A 2 niAi||c/( KlK2fc ' <T1,J2fe )||n2A2niAi) x 
K(nil^n 2 l 2 N2 uiLiS 1 QiM Ql u 2 L 2 S 2 Q 2 M Q2 LSM L Ms\ 



x a 



AM) 



(kkti) 



a (Ai) 



x a 



.(Ax) 



(«2^2)' 



nili 1 n 2 l 2 2 a[L[S[Q[MQ 1 a 2 L 2 S 2 Q 2 MQ 2 L'S'M' L M' s ). 



(kk) 



p-p 



(55) 



are schematically written down in a following as (52) and (53), by using 
expression (4.3) from Jucys and Savukynas 1973, 
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(aiJiWll K° x 4 k2) ] (k) \W1JW2f2f) = [J,f,k] 1/2 x 



x( aiJ1 pS fcl) ||aU)(a 2 i 2 ||4 K ^||a^)<; j[ j> 2 f 

k\ k 2 k 



(ta)i 



(56) 



and in the fourth case, after reversing the order of shells and altering the 
coupling of their momenta for bra and ket functions we obtain: 



{n l l^ 1 n 2 l2 2 a l L 1 S 1 Q l M Ql a 2 L 2 S 2 Q 2 M Q2 LS\\ 



,.(Ai) 



\\n 1 l^n 2 l^a[L[S[Q[M Ql a' 2 L' 2 S 2 Q' 2 M Q2 L'S' 



x a 



„(Ai) 



(kicti) 



,(A 2 ) 



x a 



.(A 2 ) 



{K2U2) 



(kk) 



[k] [L,S,L',S'\ 



1/2 



L\ L 2 L 
L[ L' 2 L' 

K\ K 2 k 



> < 



x^i/f 1 aiQiLiS-iH 
x(n 2 l 2 2 a 2 Q 2 L 2 S 2 \ 



x a 



^(Ai) 



(«ifi) 



5i S 2 5 

5; s 2 5' 

cti cr 2 A; 

7JV1 



> x 



,.(A 2 ) 



x a 



^2) 



(K2CT2) 



HmZf 1 aiQiLiSi)x 
\n 2 l^ a' 2 Q' 2 L' 2 S' 2 ), (57) 



(niZf 1 n2Z^ 2 aiLi5iQiMQ 1 a2^5 2 g2MQ 2 L 1 S , || 



,.(A 2 ) 



\\n 1 l^n2l^a' l L' l S' l Q' l M Ql a' 2 L' 2 S l 2 Q' 2 M Q2 L'S') 



x a 



-(A 2 ) 



(/tlCTl) 



x 



,.(Ai) 



x a 



^(Ai) 



(k 2 o- 2 ) 



7JV2 „,/ r ' c' /O' 



(fcfc) 



nV2 



= [fc] [L,S, 

x(n 2 l 2 2 a 2 Q 2 L 2 S 2 \\ 
ximl? 1 aiQiLiSi| 



L 2 L\ L 
L 2 L'i L' 

K\ K 2 k 



,(A2) 



x a 



(kicti) 



5" 2 S 

S' 2 S[ S' 

o\ a 2 k 

iN 2 



> x 



a (Ai) 



x a 



^(Ai) 



"I (k 2 ct 2 ) 



\\n 2 l^ a' 2 Q' 2 L' 2 S' 2 )x 
Hm/f 1 aiQiLi^i), (58) 



Prom this we conclude that in the third and fourth cases the usage of the 
tensorial expressions of only two-particle operator (8) allows us to success- 
fully exploit all the advantages of Racah algebra and quasispin formalism in 
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calculating the spin-angular parts of any two-particle operator matrix ele- 
ment. This, to our mind, not only simplifies the calculations considerably, 
by allowing to use the tables of irreducible tensors that are independent of 
shell occupation numbers, but also allows one to establish the zero matrix 
elements without performing explicit calculation. 

Meanwhile the situation is different when the last two cases are consid- 
ered, or the matrix elements between more complex configurations are to 
be established. This is related to the fact that using first (7) or second (8) 
tensorial forms the spin-angular part of matrix elements for those cases do 
not have shape of any expression (50), (51), (52) and (53). 

In the next paper we shall present a methodology that allows one to use 
efficiently the Racah algebra and quasispin formalism in a general case, too. 

6 Conclusion 

Preliminary usage of the generalized graphical method, irreducible tensorial 
form of the second quantization operators as well as of quasispin technique, 
while calculating the spin-angular parts of matrix elements of the energy 
operator, has demonstrated high efficiency to obtain in a uniform way the 
general expressions for the operators of physical quantities as well as for 
their matrix elements, covering the both cases of diagonal and non-diagonal 
ones with respect to quantum numbers of electronic configurations. There- 
fore it is fairly promising to formulate this methodology in a complete and 
consistent way for an arbitrary number of electronic shells with its successive 
implementation in the universal computer codes. 

Acknowledgements 

The authors are grateful to Professor Charlotte Froese Fischer for encourag- 
ing and valuable remarks. This work is part of co-operative research project 
funded by National Science Foundation under grant No. PHY-9501830 and 
by EURONET PECAM associated contract ERBCIPDCT 940025. 



24 



References 



Cowan R D 1981 The Theory of Atomic Structure and Spectra 
(Berkeley, CA: University of California Press) 

Froese Fischer C 1977 The Hartree-Fock Method for Atoms (New 
York: Wiley) 

Gaigalas G A 1985 Spectroscopy of autoionized states of atoms 
and ions (Moscow: Scientific council of spectroscopy) 43 (in Rus- 
sian) 

Gaigalas G A, Kaniauskas J M and Rudzikas Z B 1985 Liet. Fiz. 
Rink. (Sov. Phys. Collection) 25 3 

Gaigalas G A and Merkelis G V 1987 Acta Phys. Hungarica 61, 
111 

Jucys A P and Bandzaitis A A 1977 Theory of Angular Momen- 
tum in Quantum Mechanics (Mokslas: Vilnius) (in Russian). 

Jucys A P and Savukynas A J 1973 Mathematical Foundations 
of the Atomic Theory (Vilnius: Mokslas) (in Russian) 

Judd B R 1967 Second Quantization and Atomic Spectroscopy 
(Baltimore: John Hopkins Press) 

Lindgrcn I and Morrison M 1982 Atomic Many-Body Theory 
(Springer Series in Chemical Physics, vol. 13) (Berlin: Springer- 
Verlag) (2nd edition) 

Merkelis G V, Gaigalas G A, Kaniauskas J M and Rudzikas Z B 
1986 Izv. AN SSSR Ser. Fiz. 50 1403 

Merkelis G V, Gaigalas G A and Rudzikas Z B 1986 Liet. fiz. 
rink. (Sov. Phys. Collection) 25, 14 

Nikitin A A and Rudzikas Z B 1983 Foundations of the Theory 
of the Spectra of Atoms and Ions (Moscow: Nauka) (in Russian) 

Rudzikas Z B 1991 Comments At. Mol. Phys. 26 269 

Rudzikas Z B 1996 Theoretical Atomic Spectroscopy (Many-Electron 
Atom) (Cambridge: Cambridge University Press) (in press) 

Rudzikas Z B and Kaniauskas J M 1984 Quasispin and Isospin 
in the Theory of Atom (Vilnius: Mokslas) (in Russian) 



25 



